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ABSTRACT 

We  have  perfonned  quantum  mechanical  (QM)  dynamics  calculations  within  the 
independent-state  approximation  with  new  benchmark  triplet  A"  and  A’  surfaces  [B. 
Ramachandran  et  al.  J.  Chem.  Phys.  119,  9590  (2003).]  for  the  ro-vibronic  state-to-state 
measurements  of  the  reaction  0(3P)  +  HCl(v=2,j=l,6,9)  — »  OH(v’j’)  +  C1(2P)  [Zhang  et  al.  J. 
Chem.  Phys.  94,  2704  (1991)].  The  QM  and  experimental  rotational  distributions  peak  at  similar 
OH(j’)  levels,  but  the  QM  distributions  are  significantly  narrower  than  the  measurements  and 
previous  quasi-classical  dynamics  studies.  The  OH(low  j’)  populations  observed  in  the 
measurements  are  nearly  absent  in  the  QM  results.  We  have  also  performed  quasi-classical 
trajectory  with  histogram  binning  (QCT-HB)  calculations  on  these  same  benchmark  surfaces. 
The  QCT-HB  rotational  distributions,  which  are  qualitatively  consistent  with  measurements  and 
classical  dynamics  studies  using  other  surfaces,  are  much  broader  than  the  QM  results. 
Application  of  a  Gaussian  binning  correction  (QCT-GB)  dramatically  narrows  and  shifts  the 

1 

Approved  for  public  release;  distribution  unlimited. 


Preprint  JULY  25,  2004 


QCT-HB  rotational  distributions  to  be  in  very  good  agreement  with  the  QM  results.  The  large 
QCT-GB  correction  stems  from  the  special  shape  of  the  joint  distribution  of  the  classical 
rotational/vibrational  action  of  OH  products.  We  have  also  performed  QM  and  QCT  calculations 
for  the  transition,  O  +  HCl(v=0,T=300  K)  — >  OH(v’j’)  +  Cl  from  threshold  to  —130  kcal  mol1 
collision  energy  as  a  guide  for  possible  future  hyperthermal  O-atom  measurements.  We  find  in 
general  a  mixed  energy  release  into  translation  and  rotation  consistent  with  a  late  barrier  to 
reaction.  Angular  distributions  at  high  collision  energy  are  forward  peaked,  consistent  with  a 
stripping  mechanism.  Direct  collisional  excitation  channel  cross  sections,  O  +  HCl(v=0,T=300) 
— >  O  +  HCl(v’=l),  in  the  same  energy  range  are  large,  comparable  in  magnitude  to  the  reactive 
channel  cross  sections.  Although  the  A”  state  dominates  most  collision  processes,  above  ~48 
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kcal  mol"  ,  the  A’  state  plays  the  major  role  in  collisional  excitation. 
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I.  INTRODUCTION 

The  reaction,  0(3P)  +  HCICX’S)  — >  OH(X2n)  +  C1(2P),  is  an  important  prototype  for  heavy- 
light-heavy  (H-L-H)  systems  and  hydrogen  abstraction.  The  reaction  is  nearly  thenno-neutral, 
with  a  AH  of  -1.0  kcal  mol1  and  a  barrier  of  -10.6  kcal  mol'1  for  the  lowest  adiabatic  electronic 
state  involved.1  Experimental  rate  constant  data  are  available  from  a  number  of  initial  vibrational 
levels  of  HC1  over  a  broad  temperature  range.-'  Measurements  of  the  relative  populations  for  ro- 
vibronic  transitions  0(3P)  +  HCl(v=2,j=l,6,9)  — »  OH(v’j’)  +  C1(2P)  have  been  made  in  what  was 
the  first  fully  ro-vibronic  state-to-state  experiment.9'10  With  HCl(v=2)  levels  initially  populated, 
relatively  slow  collisions  of  O  +  HC1  occur  above  the  reaction  barrier,  and  OH(v’=0,l)  with  a 
large  range  of  OH(j’)  levels  are  observed.  Taken  together,  these  measurements  have  shown  many 
important  insights  including:  a  strong  HCl(v)  dependence  of  the  rate  constant,  a  vibrational- 
energy  conserving  tendency  (vibrational  adiabaticity),  insensitivity  of  product  rotational  level  to 
reagent  rotational  level,  a  fairly  even  partitioning  of  energy  among  translation,  vibration,  and 
rotation  of  products,  and  apparently  two  dynamical  mechanisms  that  give  rise  to  distinct 
rotational  populations. 

The  wealth  of  experimental  data  and  dynamical  information  make  this  reaction  a  test-bed  for 
dynamics  methods.  '  Using  recent,  highly  accurate  potential  surfaces,  '  agreement  between 
theory  and  measurements  has  been  very  good  for  the  rate  constant  and  the  ro-vibronic  state-to- 
state  measurements  in  general.13'14  However,  important  differences  remain.  Very  recently,  new 
benchmark  surfaces1  have  been  calculated  and  used  with  quantum  mechanical  dynamics  (QM)  to 
compute  the  thennal  rate  constant.15  Although  agreement  is  excellent  up  to  2,000  K,  large 
differences  between  these  calculations  and  measurements  at  higher  temperatures  point  to  a  need 
for  further  investigations. 

The  main  purpose  of  this  paper  is  to  report  time-independent  QM  calculations  with  the  new 
benchmark  triplet  surfaces  of  Ref.  1  for  the  ro-vibronic  state-to-state  measurements,  0(  P)  + 
HCl(v=2,j=l,6,9)  — »  OH(v’j’)  +  C1(2P),  of  Ref.  10.  The  present  calculations  are  a  stringent  test 
of  the  quality  of  the  surfaces  and  will  show  the  level  of  agreement  that  can  be  obtained  between 
highly  detailed  measurements  and  nearly  exact  theory.  It  will  also  be  the  first  time  to  our 
knowledge  that  these  measurements  have  been  modeled  with  quantum  dynamics.  A  previous 
quasi-classical  trajectory  (QCT)  study  using  the  surface  of  Ref.  11  predicted  the  peak  location 
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and  spread  of  OH(j’)  rotational  distributions  and  the  OH(v=l)/OH(v=0)  branching  ratio 
generally  very  well.11’13  For  the  surface  Ref.  12,  a  previous  QCT  study  showed  that  the  spread  of 
the  OH(j’)  rotational  distributions  were  reproduced  well,  but  the  OH(j’)  peak  locations  were  too 
high  in  energy  resulting  in  too  hot  a  rotational  distribution.14  An  important  result  from  the 
present  study  is  that  the  computed  QM  results  have  a  significantly  narrower  distribution  of 
rotational  states  than  the  measurements  and  previous  QCT  studies,11’13’14  with  the  OH(low  j’) 
populations  observed  in  the  measurements  being  nearly  absent  in  the  QM  results.  Peak  locations 
of  the  QM  OH(j’)  populations  are  quite  similar  to  the  measurements,  but  the  maxima  of  the 
relative  rotational  populations  are  much  larger  in  general  than  the  measurements  and  previous 
QCT  studies.  Possible  implications  of  the  present  QM  results  are  that  the  new  benchmark 
surfaces  are  not  adequate  to  describe  the  O  +  HCl(v=2,j)  — »  OH(v’j’)  +  Cl  transitions  and/or  the 
measurements  themselves  require  re-assessment.  It  could  also  be  that  the  QM  approach  requires 
inclusion  of  spin-orbit  coupling  effects  with  nearby  singlet  states  or  other  kinds  of  electronic 
state  coupling  not  taken  into  account  in  the  present  study. 

We  also  performed  quasi-classical  trajectory  calculations  with  the  usual  histogram  binning 
(QCT-HB)  on  these  same  benchmark  surfaces  to  evaluate  how  well  this  widely  used  method 
compares  to  exact  QM  methods  for  this  fairly  quantum  system.  The  present  QCT  calculations  on 
the  benchmark  surfaces  are  also  a  common  denominator  for  comparisons  with  QCT-HB  studies 
using  other  surfaces  so  that  differences  between  surfaces  can  be  assessed.  A  significant  finding 
of  the  present  study  is  that  the  QCT-HB  rotational  distributions  are  much  broader  than  the  QM 
results  with  the  same  benchmark  surfaces.  The  QCT-HB  results  are  similar  to  previous  QCT-HB 
studies  using  other  surfaces11’13’14  and  in  fair  agreement  with  measurements.  The  similarity  of  the 
QCT-HB  results  and  previous  classical  studies  implies  that  the  new  benchmark  potential  energy 
surfaces  are  perhaps  dynamically  similar  to  those  used  in  previous  studies,11'14  although  the 
barrier  heights  a  different. 

QCT  and  QM  calculations  for  the  same  surfaces  are  also  an  opportunity  to  examine 
techniques  that  may  correct  classical  methods.  A  simple  proposed  correction  method  is  the 
Gaussian  binning  (QCT-GB)  technique  ’  that  in  spirit  attempts  to  approximate  a  fully  semi- 
classical  approach.  Instead  of  using  all  trajectories  in  a  classical  vibrational  energy  bin  as  in 
histogram  binning  (QCT-HB),  the  QCT-GB  method  heavily  weights  trajectories  with  energies 
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near  the  exact  quantum  vibrational  level.  This  is  intuitively  appealing  and  automatically  de- 
emphasizes  trajectories  not  energetically  allowed.  It  helps  solve  the  common  problem  of 
rotational  distributions  that  are  too  broad,  extending  too  high  in  energy.  In  the  present  study  the 
QCT-GB  correction  dramatically  narrows  and  shifts  the  QCT-HB  rotational  distributions  to  an 
extent  not  before  seen  in  other  chemical  systems,  so  that  the  QCT-GB  and  QM  results  are  in  very 
good  agreement  for  all  transitions.  The  striking  consistency  of  the  QM  and  QCT-GB  results  and 
their  differences  with  measurements  highlight  the  fact  that  detailed  understanding  of  this 
chemical  system  is  still  far  away. 

As  the  effects  of  the  QCT-GB  correction  are  much  larger  and  much  different  than  in  previous 
studies,  we  investigate  the  cause  in  a  preliminary  way  by  examining  the  classical 
rotational/vibrational  action  joint  distribution  of  the  OH  products.  Future  work  will  investigate  in 
detail  the  underlying  fundamental  reasons  for  the  correction’s  apparent  success  in  O  +  HC1  and 
how  the  Gaussian-binning  correction  behaves  in  other  systems.  As  many  important  chemical 
systems  remain  out  of  reach  in  the  immediate  future  for  an  exact  QM  dynamics  treatment,  a 
reliable  correction  to  a  classical  trajectory  approach  will  be  extremely  useful. 

We  have  also  perfonned  QM  and  QCT  calculations  for  the  transition,  O  +  HCl(v=0,T=300 
K)  — >  OH(v’,j’)  +  Cl  from  threshold  to  -130  kcal  mol'1  collision  energy.  These  calculations  are 
intended  to  guide  future  hyperthermal  O-atom  measurements  centered  near  85.0  kcal  mol1  (8  km 
s'1).  Recently,  hyperthermal  O-atom  beams  have  investigated  novel  chemistries  for  a  number  of 
systems.  The  high  translational  energies  in  such  a  measurement  would  surmount  the  reaction 
barrier  for  O  +  HC1,  without  the  need  for  internal  excitation  of  the  reagent  HC1,  and  deposit  large 
amounts  of  energy  in  products.  We  report  on  the  excitation  function,  vibrational  and  rotational 
distributions  of  OH  products,  and  the  vibrationally  resolved  angular  distributions.  We  note  that 
above  -38.0  kcal  mol'1,  the  channel  O  +  HCl(v=0,T=300)  — >  OC1  +  H  is  energetically  allowed. 
The  present  surfaces  do  not  take  the  OC1  +  H  channel  into  account,  and  they  are  not  intended  to 
be  quantitative  above  -40  kcal  mol'1  even  for  the  OH  +  Cl  channel.  At  high  energies  therefore, 
the  dynamical  results  presented  here  will  only  be  qualitative.  However,  lacking  other  sources  of 
data  the  present  results  at  high  energies  will  be  an  important  starting  point  for  hyperthermal 
measurement  planning  and  for  future  theoretical  studies  in  this  energy  regime.  We  find  in 
general  a  mixed  energy  release  into  translation  and  rotation  consistent  with  a  late  barrier  to 


5 


Preprint  JULY  25,  2004 

reaction.  Above  ~35  kcal  mol-1  the  fraction  of  rotational  energy  in  OH  products  exceeds 
vibration.  The  rotational  distributions  are  extremely  hot  and  should  give  rise  to  distinct  OH  band- 
head  spectra.  Angular  distributions  at  high  collision  energy  are  forward  peaked  consistent  with  a 
stripping  mechanism.  We  have  also  performed  QCT  calculations  on  the  direct  collisional 
excitation  channel,  O  +  HCl(v=0,T=300  K)  — >  O  +  HCl(v’=l),  with  the  same  benchmark 
surfaces.  Vibrational  excitation  cross  sections  are  large,  comparable  in  magnitude  to  the  reactive 
channel  cross  sections. 

The  paper  proceeds  as  follows.  In  section  II  we  review  the  surfaces,  and  the  QM  and  QCT 
dynamical  methods  used  in  this  study.  In  section  III,  we  present  theoretical  results  for  the  0(  P) 
+  HCl(v=0,T)  — >  OH  +  Cl  rate  constant,  the  ro-vibronic  state-specific  transitions  0(  P)  + 
HCl(v=2,j=l,6,9)  ->  OH(v’j’)  +  C1(2P),  the  transitions  O  +  HCl(v=0,T=300  K)  ->  OH(v’j’)  +  Cl, 
and  the  transitions  O  +  HCl(v=0,T=300  K)  — >  O  +  HCl(v’=l),  all  with  comparisons  to 
measurements  and  other  theoretical  results  where  available.  In  section  IV  we  summarize  the 
results. 


II.  METHODS 

A.  Potential  surfaces 

3  1 

The  present  study  uses  the  two  lowest  triplet  adiabatic  surfaces  for  0(  P)  +  HC1(X  Z), 
A”  and  A’,  as  reported  in  Ref.  1  and  generated  from  a  computer  program  kindly  provided  by 
the  authors  of  that  study.  The  surfaces  are  fits  of  complete  basis  set  extrapolated  multi-reference 
configuration  interaction  calculations  with  a  large  basis  set,  and  they  are  the  most  accurate 
surfaces  available  to  date.  The  lower  A”  state  has  a  bent  (136.0  )  transition  state  geometry  with 
a  barrier  of  10.60  kcal  mol'  ,  and  the  A’  transition  state  is  linear  with  a  barrier  of  13.77  kcal 
mol'1 .  Further  details  of  the  surfaces  can  be  found  in  Ref.  1. 

All  scattering  calculations  were  performed  within  the  independent,  non-interacting  surface 
approximation.  The  repulsive  2  A’  ’  state  has  not  been  included,  and  coupling  to  singlet  surfaces 
is  ignored.  Dynamical  results  from  the  A”  and  A’  surfaces  are  weighted  by  the  following 
temperature  dependent  expression  that  attempts  to  account  for  the  spin-orbit  splitting  of  the 
overall  triplet  reagents: 
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f(T)  =  3/(5  +  3W228/r +W326/r)-  (1) 

At  300  K,  f(T)  =  0.445,  and  at  high  temperature  where  all  spin-orbit  states  become  equally 
populated  f(T)  =  1/3.  We  note  that  use  of  Eq.  (1)  for  inclusion  of  spin-orbit  degeneracy  is 
approximate  as  we  are  not  accounting  for  fine-structure  explicitly,  and  for  low  temperatures 
results  in  weights  totaling  greater  than  1.0  if  summed  over  all  three  electronic  states. 


B.  Quantum  dynamics 

28 

The  QM  calculations  were  done  using  the  quantum  reactive  scattering  code  ‘ABC’." 
Detailed  state-to-state  reaction  probabilities  P^Kv,  .,K,(E)  were  calculated  as  a  function  of  the 

total  energy  E  and  the  total  angular  momentum  J.  From  these  probabilities  the  reaction  cross 
sections  summed  over  J,  K  and  K'  for  a  particular  vj  — »  v’j’  transition  were  obtained  from  the 
standard  expression, 


°viyr(E)  = 


n 


(2  j  + 1  )kf 


■ZZ(2J  +  1)ZPiv'/U£), 


(2) 


K  J 


where  ki  is  the  wave  vector.  The  details  of  the  scattering  calculations  for  the  0(3P)+HC1  reaction 

33  •  •  15 

on  the  A"  and  A’  surfaces  are  essentially  the  same  as  in  a  recent  paper. 


For  the  O  +  HC1  (v= 2,  j  =1,6,9)  — >  OH  (v’=0,l,j’)  +  Cl  transition  calculations,  converged 
results  were  obtained  with  maximum  values  of  J  and  K(K' )  of  150  and  9,  respectively.  For  the 
coupled  channel  aspect  of  the  calculation  we  retained  all  reactant  and  product  channels  with 
internal  energies  less  than  Emax  =  55.35  kcal  mol'1  (relative  to  O  +  HC1  (rc))  and  with  diatomic 
rotational  quantum  numbers  jmax<  24.  A  total  of  2275  close-coupled  equations  were  integrated 
using  120  hyperradii  sectors  between  pmin  =1.8  a.u.  for  the  A”  surface,  1.9  a.u.  for  the  A  ’ 
surface  and  pmax  =  16.0  a.u.  for  both  surfaces.  In  order  to  simulate  the  experimentally 
measured10  product  OH  relative  v’j’  populations,  we  calculated  the  relative  populations 
following  the  approach  given  in  previous  QCT  calculations.11  Thus,  the  relative  population 
P(v'J)  is  given  by 


P(v\f)  =  F(v\f)/Ftc 


(3) 


where 
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F(v',f)  =  £  )£(2J  + 1  )Zpi, »'/*■(£«,«),  w 

i  J  K,K' 

and 

F,«=Z[m/)+F(i,/)].  e) 

/ 

In  Eq.  (4)  above,  j[Ecom )  is  the  fitted  experimental  initial  translational  energy  distribution 
function  described  below,  and  Econ  =E-Evj. 

For  the  O  +  HCl(v=0,T=300  K)  — >  OH(v’,j’)  +  Cl  transition  calculations,  we  applied  the 
standard  thermal  averaging  to  the  distribution  of  HC1  internal  states.  Due  to  the  low  temperature 
of  the  system,  the  range  of  j  and  f  is  smaller  than  for  the  O  +  HCl(v=2,  j=l  ,6,9)  calculations,  and 
the  cross  sections  were  well  converged  with  Kmax  equal  to  6.  However,  since  we  are  interested  in 
higher  collision  energies,  convergence  with  respect  to  J  for  Ecou  =  33.2  kcal  mol’1  (5  km  s’1)  was 
not  reached  until  Jmax  =  220.  Higher  collision  energies  require  too  much  computational  effort,  so 
we  stopped  at  Ecou  =33.2  kcal  mol1.  All  other  parameters  in  the  quantum  reactive  scattering 
calculations  were  the  same  as  the  calculations  on  the  O  +  HC1  (v=2,  j  =1,6,9)  — »  OH  (v’=0,l,  j’) 
+  Cl  transitions,  except  for  the  number  of  coupled  channel  equations  which  was  reduced  to  1729 
as  a  result  of  the  smaller  Kmax. 

C.  Classical  dynamics 

For  the  QCT  calculations,  we  use  standard  Monte  Carlo  techniques.-  Table  I  summarizes  the 
convergence  parameters  used  in  the  various  QCT  calculations,  where  bmax  is  the  maximum 
impact  parameter  and  Atmin  is  the  minimum  integration  time  step  for  the  variable  time  step 
integrator.  The  same  procedure  has  been  followed  to  model  the  experiments  of  Ref.  10  for  the 
transitions  the  O  +  HCl(v=2,  j  =1,6,9)  — »  OH(v’=0,l,j’)  +  Cl.  The  relative  populations  P(v’j') 
were  obtained  analogous  to  the  QM  ones,  and  follow  directly  the  approach  described  in  previous 
QCT  calculations.11  The  only  difference  from  the  QM  calculations,  is  that  instead  of  peforming 
the  collision  energy  average,  the  relative  OH(v’j’)  populations  are  obtained  by  directly  sampling 
from  the  initial  translational  energy  distribution  described  below  for  each  trajectory. 
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Table  I.  Convergence  parameters  for  QCT  calculations. 


bmax 

(a.u.) 

Atmin 

(a.u.) 

Number  of 
Trajectories 

Statistical 

Uncertainty 

0(3P)  +  HC1  — >  OH  +  Cl 
rate  constant 

4.0  ->  4.8 

1000  K  — »  2500  K 

0.01 

5.e4  — »  2.e4 
1000  K— >  2500  K 

<1  % 

0(3P)  +  HCl(v=2,j=l,6,9)  — > 
OH(v’j’)  +  C1(2P) 

5.2 

0.01 

~8.e5  for 
each  HCl(j) 

<1%  QCT-HB 
<10%  QCT-GB 

0(3P)  +  HCl(v=0,T=300)  -> 
OH(v’j’)  +  C1(2P) 

3. 4-4. 8 

0.01 

2.e5  v>  6  km  s'1 
6.e4  v<  5  km  s'1 
per  coll,  energy 

1  %(v)- 1 0%(j) 

0(3P)  +  HCl(v=0,T=300)  -> 
0(3P)  +  HCl(v’) 

3. 4-4. 8 

0.01 

2.e5  v>  6  km  s'1 
6.e4  v<  5  km  s'1 
per  coll,  energy 

l%(v) 

To  assign  final  internal  energy  states  in  the  QCT  calculations,  we  define  the  continuous 
rotational  diatomic  “quantum  number”,  j\ 

(6) 


where  Jr  is  the  magnitude  of  the  product  diatomic  angular  momentum  vector,  The  continuous 
diatomic  vibrational  quantum  number,  v  ’,  is  assigned  with  the  following  standard  semi-classical 
expression, 


7th 


VD{r) 


Jr 

2  jur2 


2’ 


(7) 


where  r±  are  the  classical  turning  points,  sint  is  the  internal  energy,  VD(r)  is  the  diatomic 

potential,  and  //  is  the  reduced  mass.  For  the  QCT-HB  calculations,  v’  and  j  ’  for  a  given 
trajectory  are  assigned  to  the  nearest  integer  bin  with  equal  weighting  wherever  they  fall  within 
the  bin.  When  these  “quantum  numbers”  are  translated  back  into  integer  bins,  they  can 
correspond  to  energetically  closed  channels,  a  well-known  problem  of  QCT-HB.  As  pointed  out 
in  Refs.  25  and  26,  one  can  lessen  this  problem  by  applying  a  Gaussian  shaped  weighting 
function  to  the  classical  trajectories  (QCT-GB)  such  that  the  trajectories  with  vibrational 
quantum  numbers  closest  to  an  integer  value  are  most  heavily  weighted,  and  those  near  the  edges 


9 


Preprint  JULY  25,  2004 


of  a  bin  are  effectively  not  considered.  The  weighting  function  used  in  the  present  work  has  the 
form. 


G(v)  = 


exp  (-(/?(v'-«))2), 


(8) 


where  n  is  the  bin  integer  value  n  =  (0,1,2,  ...)  and  ft  is  the  bin-width  parameter.  For  the  QCT- 
GB  results,  we  used  a  value  of  /?=  16.651  which  corresponds  to  a  full-width-half-maximum 
(FWHM)  of  0.1.  This  width  is  quite  narrow  and  effectively  removes  a  large  fraction  of  the 
trajectories.  The  results  were  insensitive  to  the  exact  value  of  /?,  varying  by  2-4  %  over  a  value 
corresponding  to  FWHM  from  0.05  to  0.2.  The  QCT-GB  correction  has  been  used  in  a  number 
of  recent  studies,  and  the  major  effect  has  been  to  remove  the  hotter,  energetically  closed 
rotational  populations.  We  note  that  in  spirit  the  correction  attempts  to  achieve  a  semi-classical 
approach.  However,  since  only  a  fraction  of  the  trajectories  falling  in  a  bin  contribute 
significantly,  many  more  trajectories  need  to  be  computed  for  statistical  accuracy  comparable  to 
QCT-HB.  As  an  indication  of  the  comparative  statistics  between  QCT-HB  and  QCT-GB,  Table 
II  gives  the  number  of  trajectories  for  all  reactive  transitions,  O  +  HCl(v=2,j=l,6,9)  ^ 
OH(v’=l,j’)  +  Cl,  for  the  QCT-HB  and  QCT-GB  trajectory  calculations.  Of  the  800,000  classical 
trajectories  run  for  these  transitions,  61614  contributed  to  OH  product.  Of  these  41371 
contributed  to  the  QCT-HB  result  for  v’=l  (0.5  <v’  <  1.5),  with  6163  in  the  range  0.95<v’  <  1.05 
most  important  for  QCT-GB  with  a  full-width  half-maximum  of  0.1 

Table  II.  Number  of  QCT  reactive  trajectories  for  the  transition  O  +  HCl(v=2,j=l,6,9)  ^ 
OH(v’=l,j’)  as  a  function  of  the  bin  width. 


j 

0.5<  v’  <1.5 

0.9<v’  <1.1 

0.95<v’  <  1.05 

1 

41371 

12527 

6163 

6 

10982 

2488 

1301 

9 

6898 

2349 

1145 

D.  Initial  translational  energy  distributions  for  O  +  HCl(v=2,j=l,6,9)  — >  OH(v’,j’)  +  Cl 

Modeling  of  the  state  specific  ro-vibronic  experiments  of  Ref.  10  for  the  ro-vibronic 
measurements  of  the  transitions,  O  +  HCl(v=2,j=l,6,9)  — »  OH(v’,j’)  +  Cl,  is  complicated  by  the 
fact  that  collisions  take  place  over  a  range  of  initial  translational  energies.  Ultimately,  the  center 
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of  mass  frame  translational  energy  distribution  of  the  O  +  HC1  reagents  must  account  for  the 
relative  thermal  motion  of  the  HC1  with  respect  to  the  initial  O  atom  velocities  imparted  from 
photodissociation  of  NCF.  ’  The  general  form  of  the  translational  energy  distribution  for  such 
collisions  has  been  derived  in  the  literature.  For  O  +  HC1  Ref.  14  reports,  without  supporting 
details,  a  reagent  translational  energy  distribution  that  is  Boltzmann-like  centered  at  -2.3  kcal 
mol'1  and  extending  to  -8.1  kcal  mol'1.  Since  there  remains  some  ambiguity  about  the  actual 
initial  translational  energy  distribution  for  O  +  HC1,  in  Appendix  A  we  explicitly  give  the 
functional  form  for  the  translational  energy  distribution  specifically  applied  to  the  O  +  HC1 
measurements  of  Ref.  10,  and  the  details  of  its  derivation  are  discussed.  We  present  a  Monte 
Carlo  sampling  method  to  determine  initial  collision  velocities  based  on  this  functional  form  that 
is  particularly  convenient  and  efficient  for  the  present  QCT  calculations.  The  initial  translational 
energy  distributions  obtained  are  similar  to  those  reported  in  Ref.  14. 


III.  RESULTS 

A.  Rate  constant  for  O  +  HCl(v=0)  — >  OH  +  Cl 

Figure  1  shows  the  present  calculated  QCT-HB  and  QM  results15  along  with  recent 
experimental  results  between  1,000  K  and  -3,000  K.  Both  QCT-HB  and  QM  results  include  the 
A”  and  A’  surfaces.  As  shown  previously,  the  QM  results  are  in  excellent  agreement  with 
measurements  up  to  -2,000  K.  Above  this  temperature,  it  was  shown  that  the  OC1  channel  is 
likely  not  contributing,  and  the  differences  between  theory  and  measurement  remain  an  open 
issue.  The  QCT-HB  results  are  about  a  factor  of  two  lower  than  the  QM  results  at  1,000  K,  but 
by  2,500  K  they  are  nearly  identical.  This  is  consistent  with  a  previous  study  using  other  surfaces 
at  lower  temperatures14  and  with  the  expectation  that  tunneling  is  becoming  less  important  with 
increasing  energy.  These  calculations  help  confirm  that  the  O+HCl  potentials  are  being  used 
consistently  in  the  quantum  and  classical  scattering  calculations. 

B.  Cross  sections  and  relative  populations  for  O  +  HCl(v=2,j=l,6,9)  — >  OH(v’j’)  +  Cl 

We  wish  to  model  the  measurements  of  Ref.  10  for  the  transitions,  O  +  HCl(v=2,j=l,6,9)  — > 
OH(v’j’)  +  Cl.  As  discussed  earlier,  the  initial  collision  energies  have  a  Boltzmann-like  spread 
from  0  up  to  -8  kcal  mol"1.  To  understand  the  impact  of  this  energy  distribution  and  for  a  more 
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straightforward  comparison  between  QM  and  QCT  results,  in  Figure  2  we  show  the  energy 
dependence  of  the  cross  sections  for  the  present  calculated  QCT-HB  and  QM  results  for  the 
transitions  O  +  HCl(v=2,j=l,6,9)  — »  OH(v’=0,l)  +  Cl.  For  comparison  purposes  these  results  use 
only  the  A”  surface.  (The  A’  surface  makes  a  relatively  small  contribution  to  the  cross 
section.)  Except  very  near  threshold,  the  agreement  between  QCT-HB  and  QM  results  is 
remarkably  good  including  the  local  structure  between  0.1  and  0.2  eV  where  QCT  calculations 
were  done  for  HCl(v=2,j=6,9).  (We  wish  to  emphasize  that  the  QCT  cross  sections  at  fixed 
energy  presented  in  Figure  2  were  not  used  to  generate  results  for  comparison  purposes  with  the 
measurements  of  Ref.  10.  To  compute  relative  populations,  QCT  results  were  obtained  by 
directly  sampling  from  the  initial  translational  energy  distribution  discussed  earlier.)  We  note 
that  as  in  Ref.  14  Fig.  11,  the  HCl(v=2,j=l)  cross  section  is  much  larger  than  the  other 
transitions.  However,  all  cross  sections  presented  here  are  about  a  factor  of  2  larger  than  those  in 
Ref.  14.  The  relative  magnitude  ordering  of  the  HCl(v=2,j)  cross  sections  of  the  present  results 
are  much  different  than  those  using  the  surface  of  Ref.  11. 

We  note  that  the  collision  energy  distribution  in  the  measurements  of  Ref.  10  are  Boltzmann- 
like  with  a  maximum  near  0.1  eV,  but  the  calculated  energy  dependent  cross  sections  continue  to 
rise  rapidly  after  this  energy.  The  net  effect  is  that  collisions  near  0.2  eV  appear  to  be  the  most 
important.  Furthermore,  we  have  found  that  results  for  a  fixed  collision  energy  around  0.2  eV  are 
quite  similar  to  results  taking  into  account  the  full  initial  translational  energy  distribution. 

Figure  3  shows  the  present  theoretical  results,  QM,  QCT-HB  (histogram  binning),  and  QCT- 
GB  (Gaussian  binning),  along  with  the  measurements  of  Ref.  10  for  the  transitions  O  + 
HCl(v=2,j=l,6,9)  — »  OH(v’=0,l,j’)  +  Cl  in  terms  of  relative  populations  for  each  initial 
HCl(v=2,j)  state.  The  theoretical  results  use  the  full  initial  translational  energy  distribution 
discussed  earlier  and  both  A”  and  A’  states. 

A  significant  result  of  the  present  calculations  as  shown  in  Figure  3  is  that  the  QM  relative 
populations  are  much  narrower  and  larger  at  the  peak  than  the  measurements,  with  the  prominent 
OH(low  j’)  populations  in  the  measured  results  nearly  absent  in  the  QM  results.  The  peak 
locations  of  the  rotational  populations  agree  very  well  between  the  QM  and  experimental  results, 
especially  for  OH(v’=l).  The  differences  between  QM  and  experimental  results  are  consistent 
across  the  HCl(v=2,j)  initial  states  and  extend  far  beyond  measured  error  bars.  This  may  mean 
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that  the  new  benchmark  surfaces  are  not  adequate  to  describe  the  O  +  HCl(v=2,j)  — >  OH(v’j’)  + 
Cl  transitions  and/or  the  measurements  themselves  require  re-assessment.  It  may  also  point  to  the 
importance  of  spin-orbit  coupling  effects  with  nearby  singlet  states  or  other  kinds  of  electronic 
state  interactions  not  taken  into  account  in  the  present  study.  We  note  for  example  that  Ref.  10 
measures  the  OH  II3/2  component  only,  and  this  may  have  to  be  explicitly  accounted  for  in  the 
modeling. 

As  a  point  of  reference,  the  present  QCT-HB  results  are  in  fair  agreement  with  the 
measurements,  although  consistently  extending  too  high  in  OH(j’)  and  predicting  an  OH(j’) 
maximum  that  is  too  low  in  OH(j’).  The  QCT-HB  results  have  the  same  level  of  agreement  with 
earlier  QCT  results  using  a  different  surface.11’13  The  largest  difference  between  the  present 
QCT-HB  results  and  measurements  occur  for  HCl(v=2,j=6). 

Another  significant  finding  of  the  present  study  is  that  the  QCT-GB  results  are  dramatically 
different  from  the  QCT-HB  results  and  closely  resemble  the  QM  results  for  all  transitions.  The 
QCT-GB  correction,  in  addition  to  lessening  the  very  high  j’  populations  as  seen  in  earlier 
studies  on  other  systems,'  severely  diminishes  the  low  j’  populations,  narrowing  and  shifting  the 
j’  distribution  to  match  the  QM  results  very  well.  Since  the  QCT-GB  correction  seen  here  is  far 
larger  and  different  than  seen  in  other  systems,  we  have  examined  the  classical  joint 
vibrational/rotational  (continuous)  distributions  to  understand  the  underlying  mechanism  for  its 
apparent  success.  In  Figure  4  we  show  a  contour  plot  of  the  joint  probability  as  a  function  of  the 
continuous  classical  vibrational  and  rotational  quantum  numbers  using  trajectories  for  the 
transition  O  +  HCl(v=2,j=l)  — >  OH(v’j’)  +  Cl.  The  contours  have  a  distinct  ridge-like  shape 
along  a  line  of  high-jci,  low-vci  and  extend  to  low-jci,  high-vci,  cutting  off  on  the  high  energy  side 
at  the  maximum  available  product  energy.  The  ridge  extends  from  the  maximum  near  vci=1.25, 
jci=7.5  toward  the  low-vci,high-jci  corner,  and  appreciably  stretches  into  the  vci=2  bin. 

For  O  +  HCl(v=2,j=l)  — >  OH(v’=0,j’)  +  Cl  transitions  for  example,  QCT-HB  will  use  all 
trajectories  with  vci  between  -0.5  and  0.5  with  equal  weighting,  while  QCT-GB  will  weight  the 
trajectories  with  a  Gaussian  function  such  that  those  trajectories  within  the  shaded  band  around 
vci  =0.0  will  be  very  heavily  weighted.  Because  of  the  ridge-like  shape  of  the  probability  and  its 
orientation,  it  is  clear  how  application  of  the  Gaussian  weighting  narrows  the  rotational 
distribution.  Because  within  a  full  bin  the  probability  is  weighted  toward  high  vci  and  low-jci, 
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application  of  the  weighting  function  removes  a  great  deal  of  high-v(low-j)  populations.  It  is  also 
interesting  that  weighting  recovers  the  QM  relative  population  magnitudes.  Although  promising, 
this  work  is  preliminary,  and  we  do  not  report  a  Gaussian  binning  correction  to  other  results  in 
this  paper.  But  we  do  note  that  the  Gaussian  binning  apparently  has  only  a  small  affect  on 
vibrational  and  total  cross  sections.  Future  work  will  investigate  in  detail  the  underlying 
fundamental  reasons  for  the  correction’s  apparent  success  in  O  +  HC1  and  how  the  Gaussian- 
binning  correction  behaves  in  O+HCl  transitions  and  in  other  systems. 

Table  III  gives  total  and  fractional  product  energies  as  a  function  of  initial  HCl(j)  level.  The 
experiment  entries  are  taken  from  Ref.  10.  We  adopt  their  convention  for  defining  energies  for  a 
more  direct  comparison:  <Etot>  is  defined  as  the  sum  of  initial  vibrational  energy,  average  initial 
translational  energy,  initial  rotational  energy  and  the  heat  of  fonnation,  <EV>  is  the  average  final 
OH  vibrational  energy  with  respect  to  the  OH  zero  point,  <Er>>  is  the  average  final  OH  rotational 
energy,  and  the  average  final  translational  energy  is  defined  <ET’>  =  <Etot’>  -  (<EV’>  +  <Er’>). 
We  also  adopt  the  convention  for  the  average  fractional  final  energies,  <fv’,r’,T,:>  = 
<Ev’,Er’,Er>/<Etot’>.  The  table  entries  reflect  results  evident  from  the  figures.  Since  the  low 
OH(j’)  populations  observed  in  the  measurements  are  nearly  absent  in  the  QM  and  QCT-GB 
results,  the  average  rotational  energy  fraction  in  the  QM  and  QCT-GB  results  are  larger  than  the 
measurements  for  all  transitions.  In  fact,  for  all  transitions  the  QM  fractional  rotational  energies 
are  larger  than  the  QM  fractional  vibrational  energies,  which  is  opposite  to  what  is  observed 
experimentally.  The  fractional  amount  of  internal  energy  is  greater  for  all  theoretical  results 
compared  to  measurements,  and  the  partitioning  of  energy  between  vibration  and  rotation  in  the 
theoretical  results  varies  more  as  a  function  initial  HCl(j)  level  than  the  measurements.  The 
QCT-HB  average  rotational  fractions  are  in  fair  agreement  with  the  experimental  results. 
Including  contributions  from  the  OH(v’=2)  QCT  results  makes  the  fraction  of  energy  going  into 
vibration  much  larger,  as  expected  from  examination  of  the  joint  classical  vibration/rotation 
actions. 
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Table  III.  Average  total  and  fractional  energy  partitioning  for  OH(v’j’)  product  for  transitions 
HCl(v=2,j)  — >  OH(v’j’)  +  Cl.  Main  entries  include  OH(v’=0,l)  states  and  for  QCT  results  are 
obtained  by  assigning  energies  after  quantum  binning.  QCT  entries  in  parentheses  include  all 
transitions  to  v’=0,l,  and  2  OH  states,  and  they  were  obtained  by  assigning  energies  of  classical 
trajectories  directly,  before  binning  into  quantum  OH(v’j’)  final  states.  All  energies  are  in  kcal 
mol'1. 


Experiment 
Ref.  10 

QM 

QCT-HB 

QCT-GB 

<Etot’> 

18.6 

<EV> 

7.3 

8.2 

8.3 

8.2 

<Er> 

5.9 

7.3 

6.2 

7.4 

HC1G=1) 

<Ex’> 

5.4 

3.1 

4.0 

3.1 

<fv’> 

0.395 

0.439 

0.449(0.547) 

0.441(0.440) 

<fr> 

0.320 

0.495 

0.331(0.272) 

0.395(0.361) 

<fr> 

0.285 

0.166 

0.220(0.181) 

0.164(0.199) 

<Etot*> 

19.8 

<EV> 

7.6 

8.0 

9.2 

9.2 

<Er-> 

6.4 

8.5 

5.4 

7.2 

HCl(j=6) 

<Et> 

5.8 

3.3 

5.2 

3.5 

<fv’> 

0.386 

0.403 

0.464(0.622) 

0.463(0.504) 

<fj’> 

0.325 

0.431 

0.273(0.210) 

0.363(0.310) 

<fT’> 

0.289 

0.166 

0.263(0.168) 

0.174(0.186) 

<Etot*> 

21.2 

<EV> 

8.2 

8.5 

8.9 

8.7 

<Er> 

7.0 

8.6 

7.0 

8.8 

HCl(j=9) 

<Er> 

6.0 

4.1 

5.3 

3.8 

<fv-> 

0.386 

0.402 

0.421(0.611) 

0.408(0.602) 

<fr> 

0.331 

0.406 

0.330(0.212) 

0.414(0.242) 

<fT>> 

0.284 

0.192 

0.249(0.177) 

0.178(0.156) 

Table  IV  gives  the  cross  section  summed  over  OH(v’=0,l,j’)  as  a  function  of  HCl(j)  level. 
Although  the  QCT  and  QM  results  are  different  in  magnitude,  they  show  a  very  similar  decrease 
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in  cross  section  as  a  function  of  the  initial  HCl(j):  a  factor  of  3.4  for  QCT-HB  and  5.3  for  QM  in 
going  from  HCl(j=l)  to  HCl(j=9).  In  contrast,  the  experimental  results  increase  as  a  function  of 
HCl(j)  by  a  factor  of  1.5  ±0.5  in  going  from  HCl(j=l)  to  HCl(j=9).10 


Table  IV.  Total  cross  section  summed  over  all  final  states  OH(v’j’)  as  a  function  of  HCl(j)  level 
for  the  transition  HCl(v=2,j)  — »  OH(v’,j’)  +  Cl. 


j 

QCT-HB 

a(a.u.) 

QM 

a(a.u.) 

1 

1.294 

2.869 

6 

0.440 

0.685 

9 

0.386 

0.546 

The  vibrational  branching  ratios,  a(v’=l)/a(v’=0)  are  shown  in  Table  V  as  a  function  of 
HCl(j)  level.  For  HCl(j=T)  all  theoretical  results  are  in  fair  agreement  with  each  other  and  they 
are  larger  than  the  experimental  results.  The  QM  and  experimental  results  are  in  fair  agreement 
for  HCl(j=6,9).  However,  the  QCT  results  are  all  consistently  much  larger  than  the  experiment, 
especially  for  HCl(j=6),  although  the  QCT-GB  correction  does  improve  agreement  with  the  QM 
results.  The  under-estimate  of  the  v’=0  level  contribution  in  the  QCT-HB  results  is  not  fully 
recovered  in  the  QCT-GB  correction. 


Table  V.  Vibrational  Branching  Ratios  OH(v’=l)/OH(v’=0)  for  O  +  HCl(v=2,j)  — >  OH(v’,j’)  + 
Cl 


j=l 

j=6 

j=9 

Experiment 

2.6±0.1 

3.0±0.1 

4.1±0.2 

QM 

4.0 

3.6 

5.1 

QCT-HB 

4.5 

9.1 

7.0 

QCT-GB 

4.1 

8.7 

5.6 
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C.  Cross  sections  for  O  +  HCl(v=0,T=300  K)  — >  OH(v’j’)  +  Cl 

As  a  guide  to  possible  future  hyperthermal  O-atom  beam  measurements  centered  near  ~8 
km  s'1  (85.0  kcal  mol'1),  which  surmount  the  reaction  energy  barrier  through  reagent  translation 
as  opposed  to  internal  energy,  we  present  results  for  the  transition  O  +  HCl(v=0,T=300  K)  — » 
OH(v’j’)  +  Cl  from  threshold  to  10  km  s'1  collision  velocity  (-133  kcal  mol’1).  This  is  the  first 
time  to  our  knowledge  that  results  for  this  system  have  been  presented  above  ~20  kcal  mol1.  The 
upper  energy  range  of  these  calculations  is  far  beyond  the  energy  threshold  for  the  OC1  +  H 
channel  at  (-38  kcal  mol'1)  and  the  reported  valid  range  of  the  potential  surfaces  even  for  the  OH 
+  Cl  channel.  The  present  results  above  ~40  kcal  mol  1  (~5.5  km  s'1)  use  what  is  really  an 
extrapolation  of  the  fitted  surfaces,  and  so  there  is  a  large  uncertainty  associated  with  any 
dynamical  results  generated  from  them.  However,  the  present  calculations  yield  important 
qualitative  information  and  serve  as  a  baseline  for  future  work. 

Figure  5  shows  results  for  the  QM  and  QCT-HB  reactive  cross  section  for  O  + 
HCl(v=0,T=300  K)  — »  OH(v’)  +  Cl  as  a  function  of  relative  collision  velocity.  The  QM  results 
include  only  the  A”  surface  and  the  QCT-HB  results  include  both  surfaces.  Figure  5a  shows  the 

3  3 

QCT-HB  total  reactive  cross  section,  the  A”  and  A’  component  contributions,  and  the  QM 

-3 

results  for  the  A”  surface  alone.  There  is  excellent  agreement  between  the  QM  and  QCT-HB 
results.  They  show  a  steep  rise  from  threshold  and  a  leveling  off  of  the  total  cross  section  which 
is  dominated  by  the  A”  contribution.  Above  ~4  km  s'  (21.3  kcal  mol"  )  the  cross  section 
remains  constant  over  a  large  velocity  range.  Figure  5b  shows  the  vibrationally  resolved  cross 
sections  as  a  function  of  velocity.  There  is  generally  good  agreement  between  QM  and  QCT-HB 
results,  but  the  v’=0  cross  sections  are  underestimated  and  the  v’=l  are  over-estimated  by  the 
QCT-HB  method,  consistent  with  our  findings  in  Table  V.  There  is  a  large  falloff  in  the  cross 
section  with  vibrational  quanta,  especially  near  threshold,  with  most  of  the  product  energy  going 
into  translation  and  rotation.  This  is  consistent  with  a  late  barrier  picture  of  the  energy  release. 

Figure  6  shows  the  OH(j’)  cross  sections  for  OH(v’=0,l)  at  collision  velocities  of  4,  6,  and  8 
km  s'1  (21.3,  47.8,  and  85.0  kcal  mol"1,  respectively).  QCT-HB  results  are  shown  for  all 
velocities,  and  QM  results  are  shown  at  4  km  s’1.  The  rotational  distributions  peak  approximately 
between  j’=10  and  j ’=20  (~6  and  ~23  kcal  mol’1)  and  extend  up  to  j’=35  (68  kcal  mol'1)  at  8  km 
s'1.  At  4  km  s'1,  agreement  between  QM  and  QCT-HB  results  is  very  good  with  the  QM  results 
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slightly  narrower,  but  with  the  same  peak  locations.  The  QCT-HB  results  extend  farther  in 
OH(j’)  than  the  QM  results,  behavior  that  could  perhaps  be  corrected  with  Gaussian  binning. 
Unlike  the  previous  O  +  HCl(v=2,j)  QM  and  QCT-GB  results,  however,  low  OH(j’)  populations 
are  significant,  and  it  would  be  interesting  to  see  if  the  QCT-GB  correction  can  be  effectively 
applied  here.  Figure  7  shows  the  average  rotational  quantum  number  of  the  OH(j’)  product  as  a 
function  of  collision  velocity.  There  is  good  agreement  between  the  QCT-HB  and  QM  results 
except  near  threshold,  showing  the  steep  increase  in  OH(j’)  with  velocity.  Figure  8  shows  the 
fractional  energy  release  in  products  as  a  function  of  collision  velocity.  Above  5  km  s'1,  the 
fraction  of  energy  in  rotation  exceeds  vibration.  The  high  OH(j’)  populations  will  produce  a 
distinct  band-head  OH  spectrum  for  the  nascent  products,  something  that  should  be 
experimentally  observable.  The  strong  translational  and  rotational  energy  release  in  this  reaction 
is  consistent  with  the  late  barrier  in  the  A”  potential  energy  surface.  Finally,  Figure  9  shows  the 
OH(v’)  angular  distribution  as  a  function  of  the  center  of  mass  scattering  angle.  The  dominant 
v’=0  and  v’=l  cross  sections  are  quite  forward  peaked  indicating  a  stripping  mechanism. 

D.  Collisional  excitation  of  HC1:  O  +  HCl(v=0,T=300  K)  — >  O  +  HCl(v’=l) 

Hyperthennal  O  atom  observations  may  also  contain  a  substantial  collisional  excitation 
component.  To  investigate  this  possibility  we  have  performed  QCT-HB  calculations  on  the 
transition  O  +  HCl(v=0,T=300  K)  — >  O  +  HCl(v’=l).  The  results  are  shown  in  Figure  10  broken 
down  into  A”  and  A’  components.  Near  threshold,  the  A”  surface  contribution  dominates  and 
the  cross  sections  are  relatively  small.  Above  6  km  s'1,  however,  the  3A’  surface  contribution 
becomes  most  important  in  what  is  likely  a  different  excitation  mechanism.  It  could  be  that 
trajectories  that  are  reactive  on  the  A”  surface  create  vibrational  excitation  very  effectively  on 
the  A’  surface.  At  these  high  velocities,  the  total  cross  section  begins  a  steep  rise  making 
collisional  excitation  channel  comparable  in  magnitude  to  the  reactive  channel. 


IV.  SUMMARY 

We  have  performed  time-independent  quantum  mechanical  (QM)  dynamics  calculations  with 
the  new  benchmark  triplet  surfaces  of  Ref  1  for  the  ro-vibronic  state-to-state  scattering,  0(  P)  + 
HCl(v=2,j=l,6,9)  — »  OH(v’j’)  +  C1(2P),  of  Ref.  10  and  compared  to  experiment.  Although  the 
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OH(j’)  peak  locations  agree  well  between  QM  results  and  experiment,  the  QM  results  have  a 
significantly  narrower  distribution  of  rotational  states  than  measurements  and  previous  classical 
dynamics  studies.  The  OH(low  j’)  populations  observed  in  the  measurements  are  nearly  absent  in 
the  QM  results.  Possible  implications  of  the  present  QM  results  are  that  the  new  benchmark 
surfaces  are  not  adequate  to  describe  the  O  +  HCl(v=2,j)  — »  OH(v’j’)  +  Cl  transitions  and/or  the 
measurements  themselves  require  re-assessment.  It  could  also  be  that  the  QM  approach  requires 
inclusion  of  spin-orbit  coupling  effects  with  nearby  singlet  states  or  other  kinds  of  electronic 
state  interactions  not  taken  into  account  in  the  present  study. 

We  have  also  performed  quasi-classical  trajectory  with  histogram  binning  (QCT-HB) 
calculations  on  these  same  benchmark  surfaces.  The  QCT-HB  rotational  distributions,  which  are 
qualitatively  consistent  with  measurements  and  classical  dynamics  studies  using  other  surfaces, 
are  much  broader  than  the  QM  results.  Application  of  a  Gaussian  binning  correction  (QCT-GB) 
dramatically  narrows  and  shifts  the  QCT-HB  rotational  distributions  to  be  in  very  good 
agreement  with  the  QM  results.  The  cause  of  the  large  QCT-GB  correction  stems  from  the 
special  shape  of  the  classical  rotational/vibrational  action  joint  distribution  of  OH  products. 
Future  work  will  investigate  in  detail  the  underlying  fundamental  reasons  for  the  correction’s 
apparent  success  in  O  +  HC1  and  how  the  Gaussian-binning  correction  behaves  in  other  systems. 
The  striking  consistency  of  the  QM  and  QCT-GB  results  and  their  differences  with 
measurements  highlight  the  fact  that  detailed  understanding  of  this  chemical  system  is  still  far 
away. 

We  have  also  performed  QM  and  QCT  calculations  for  the  reactive  transitions,  O  + 
HCl(v=0,T=300  K)  — »  OH(v’,j’)  +  Cl  from  threshold  to  -130  kcal  mol'1  collision  energy  as  a 
guide  for  possible  future  hyperthermal  O-atom  measurements.  Above  -40  kcal  mol1  we  expect 
only  qualitative  results  as  this  is  beyond  the  reported  valid  range  of  the  surfaces.  However,  these 
results  will  be  an  important  baseline  for  future  measurements  and  theoretical  studies.  We  find  in 
general  a  mixed  energy  release  into  translation  and  rotation  consistent  with  a  late  barrier  to 
reaction.  The  rotational  distributions  are  extremely  hot  and  should  give  rise  to  distinct  OH  band- 
head  spectra.  Angular  distributions  at  high  collision  energy  are  forward  peaked  consistent  with  a 
stripping  mechanism.  We  have  also  performed  QCT  calculations  on  the  direct  collisional 
excitation  channel,  O  +  HCl(v=0,T=300  K)  O  +  HCl(v’=l),  in  the  same  energy  range  with  the 
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same  benchmark  surfaces.  Vibrational  excitation  cross  sections  are  large,  comparable  in 
magnitude  to  the  reactive  channel  cross  sections.  Although  the  A”  state  dominates  most 
collision  processes,  above  ~48  kcal  mol"  ,  the  A’  state  plays  the  major  role  in  collisional 
excitation. 
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APPENDIX  A.  O+HCl  translational  energy  distribution  for  the  measurements  of  Ref.  10. 

In  this  appendix,  we  review  the  derivation  of  the  center  of  mass  frame  reagent  translational 
energy  distribution  for  the  ro-vibronic  state-specific  measurements  of  Ref.  10,  O  + 
HCl(v=2,j=l,6,9)  — »  OH(v',j')  +  Cl.  The  reagent  O-atoms  are  initially  formed  from 
photodissociation  of  NO2.  The  resulting  measured  O  atom  velocity  distribution  of  Ref.  32 
(represented  by  the  sum  of  the  two  Gaussians)  is  shown  in  Figure  Al.  The  O  atom  velocity  in 
this  figure  is  in  the  center  of  mass  frame  of  the  NO?  precursor.  The  lower  velocity  Gaussian- 
shaped  peak  corresponds  to  formation  of  NO(v=l)  and  the  higher  velocity  peak  to  formation  of 
NO(v=0).  The  explicit  fonn  used  to  model  this  distribution  is, 

f(u)  =  Aexp(-(a(u  -iq))2)  +Bvxp(-(b(u  -u2))2) ,  (Al) 

where  iii  =  890  m  s"1,  U2  =  1400  m  s"1,  a  =  4.6353e-3  s  m"1,  b  =  4.6353e-3  s  m"1,  A  =  0.847,  and 
B=  0.997. 

To  model  the  O  +  HC1  measurements  of  Ref.  10,  we  need  the  O  +  HC1  center  of  mass 
frame  relative  velocity  distribution.  First,  the  velocity  in  the  NO?  frame  is  averaged  over  the 
thennal  distribution  of  the  NO2  precursor,  which  results  in  an  O  atom  distribution  in  the 
laboratory  frame  (Eq.  1  of  Ref.  30)  for  a  given  speed,  say  v.  Then  v  is  averaged  over  the  HC1 
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thermal  distribution,  as  originally  done  in  Ref.  33  (see,  also  Eq.  2  of  Ref.  30).  The  resultant 
conditional  relative  velocity  distribution  is  then  averaged  over  the  O  atom  distribution  in  the 
laboratory  frame  as  described  by  Eq.  1  of  Ref.  30,  resulting  in  a  relative  velocity  (speed) 
distribution  in  the  O  +  HC1  frame  for  a  single  value  of  u  (Eq.  3  of  Ref.  30).  If  the 
photodissociation  process  resulted  in  a  single  velocity  value,  we  would  simply  use  Eq.  3  Ref.  30. 
However,  since  we  have  a  distribution  of  velocities  given  by  Ref.  32.,  Eq.  3  of  Ref.  30  is  then 
averaged  over  the  measured  O  atom  distribution  shown  above.  The  resulting  expression  is  given 
by  Eq.  (8)  of  Ref.  31,  which  reduces  to  Eq.  3  of  Ref.  30  if  the  experimental  distribution  from  the 
photodissociation  process  is  given  by  5(v-vexc)- 

Explicitly,  for  the  processes,  AB  +  hv  — >  A  +  B  and  A  +  C  — >  products,  the  final  expression 
for  the  A  +  C  velocity  distribution  in  the  center  of  mass  frame  used  in  the  present  QCT 
calculations  is: 


f  (y  con )  ~ 


/  y/2 

L 

Kn  j 


J  duf  ( u)i 


coll 


{exp  {-y(u  -  vcollf  -  exp  (~y(u  +  vcoll  )2 }. 


(A2) 


m  in  . 

Her ey  = - £ — - ,  u  is  the  velocity  (measured)  in  the  AB  center  of  mass  frame,  f(u)  is  the 

ma  +  mh  +  mc 

velocity  distribution  given  in  Eq.  Al,  and  vcou  is  the  velocity  in  the  A  +  C  center  of  mass  frame. 
In  the  present  modeling,  NCfi  +  hv  — >  O  +  NO  and  O  +  HC1  — >  OH  +  Cl,  so  that  A  =  O,  B=NO, 
and  C  =  HC1.  Eq.  A2  is  particularly  convenient  for  selecting  initial  O  +  HC1  velocities  for  QCT 
calculations.  In  the  present  QCT  calculations  we  select  the  initial  O  +  HC1  velocities  in  two 
steps: 

1)  Select  u  from  duf(u)u 1 . 


2)  Given  u,  select  vcoii  from 


coll 


{exp  {-y{u  -  v 


coll 


Y  -  exp(-y(u  +  vcoll) 


The  resulting  O+HCl  relative  speed  distribution  from  Eq.  A2  is  shown  in  Figure  A2  as 
the  fo-HCi(v)  curve  (solid  line).  Also  shown  is  the  kernel  distribution,  Eq.  Al,  times  the  velocity 
squared  as  a  reference  (dashed  line).  Figure  A3  shows  the  relative  translational  energy 
distribution  for  O+HCl  center  of  mass  collisions.  To  judge  the  reasonableness  of  the 
distributions,  consider  the  approximate  FWHM  for  a  (single)  relative  collision  energy  of  ~0.1 
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eV.  Eq.  (4)  of  Ref.  30  estimates  the  FWHM  to  be  0.094  eV.  The  two  maxima  in  the  measured 
O-atom  distribution  of  Ref.  32  are  separated  by  -0.067  eV.  It  is  not  surprising  that  the  thermal 
averaging  washes  out  the  structure  in  the  measured  O  atom  distribution. 
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Figure  Captions: 

Figure  1.  Thermal  rate  constant  for  O  +  HCl(v=0)  — >  OH  +  Cl  from  1,000-3,000  K.  □, 

measurements  of  Ref  3;  x,  measurements  of  Ref  2; - ,  QM  results  of  Ref.  15,  — • — ,  present 

QCT-HB  results. 

Figure  2a-c.  Vibrationally  resolved  cross  sections  as  a  function  of  energy  for  the  transitions  O  + 
HCl(v=2,j=l,6,9)  — >  OH(v’=0,l)  +  Cl.  □,  QM  results  for  v’=0;  |,  QM  results  for  v’=l;  <>,  QCT- 
HB  results  for  v’=0;  ♦,  QCT-HB  results  for  v’=l. 

Figure  3a-c.  Relative  populations  of  OH(v’j)  as  a  function  of  OH(j’)  level  for  the  transitions  O  + 
HCl(v=2,j=l,6,9)  — »  OH(v’=0,l)  +  Cl.  °,  Experimental  results  of  Ref.  10  for  v’=0;  •, 

Experimental  results  of  Ref.  XXX  for  v’=l;  □,  QM  results  for  v’=0;  |,  QM  results  for  v’=l;  «, 
QCT-HB  results  for  v’=0;  ♦,  QCT-HB  results  for  v’=l;  *,  QCT-GB  results  for  v’=0;  *,  QCT-GB 
results  for  v’=l. 

Figure  4.  Joint  probabilities  of  the  continuous  QCT  rotational  (jd)  and  vibrational  (vci)  action  for 
the  transition  O  +  HCl(v=2,j=l)  — »  OH(v’j’)  +  Cl.  The  contours  (0.9,  0.7,  0.5,  0.3,  0.1,  and  0.05) 
are  normalized  such  that  the  maximum  is  1.0  (indicated  by  an  *).  Results  include  contributions 
from  both  the  A”  and  A’electronic  states. 


Figure  5.  Total  and  vibrationally  resolved  cross  sections  for  O  +  HCl(v=0,T=300)  — >  OH  +  Cl  as 
a  function  of  relative  collision  velocity,  (a)  Total  cross  sections  for  O  +  HCl(v=0,T=300)  — >  OH 
+  Cl  and  their  electronic  components.  ♦,  QCT-HB  total  cross  section  including  the  A”  and  A’ 
states,  a,  QM  results  for  the  3A”  state  up  to  5  km  s'1;  «,  QCT-HB  results  for  the  3A”  state; 
QCT-HB  results  for  the  A’  state;,  (b)  Total  vibrationally  resolved  cross  sections  for  O  + 
HCl(v=0,T=300)  — >  OH(v’)  +  Cl  and  their  electronic  components.  ♦,  QCT-HB  total  cross 
including  A”  and  A’  states  for  v’=0,l„  and  2;  «,  QCT-HB  total  cross  section  only  including  the 
A”  state  for  v’=0,l,  and  2;  □,  QM  total  cross  section  for  the  A”  state  only  up  to  5  km  s'  for 
v’=0,l,  and  2. 
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6a-c.  Rotationally  resolved  cross  sections  versus  OH(j’)  for  the  transitions  O  +  HCl(v=0,T=300) 
— »  OH(v’=0,l,j’)  +  Cl  at  4,  6  and  8  km  s'1,  respectively.  □,  QM  results  at  4  km  s'1  for  v’=0;  |, 
QM  results  at  4  km  s'1  for  v’=l;  «,  QCT-HB  results  at  4,  6,  and  8  km  s'1  for  v’=0;  ♦,  QCT-HB 
results  at  4,  6,  and  8  km  s'1  for  v’=l. 

7.  Average  OH(j’)  quantum  number  versus  relative  collision  velocity  for  the  transitions  O  + 
HCl(v=0,T=300)  — >  OH(v’j’)  +  Cl  summed  over  v’.  1,  QM  results;  ♦,  QCT-HB  results. 

8.  Fractional  energy  disposal  into  OH  products  for  the  QCT-HB  results  versus  relative  collision 
velocity  for  the  transition  O  +  HCl(v=0,T=300)  — >  OH(v’j’)  +  Cl.  ^  ,  translational  energy;  «, 
vibrational  energy;  ♦,  rotational  energy. 

9.  Vibrationally  resolved  angular  differential  cross  section,  s i n ' 1  ( 0 )  da/d0,  versus  center  of  mass 
angle,  0,  at  8  km  s'1  for  the  QCT-HB  results  for  the  transitions  O  +  HCl(v=0,T=300)  — >  OH(v’) 
+  C1. 

10.  Total  and  electronic  state  component  collisional  excitation  cross  sections  versus  relative 
collision  velocity  for  the  QCT-HB  results  for  the  transition  O  +  HCl(v=0,T=300)  — >  HCl(v’=l)  + 
Cl.  ♦  , total  cross  section;  <>,  A”  component;  x  A’  component. 


Appendix  A.  Figure  Captions: 

Figure  Al.  Fit  of  the  O  atom  speed  distribution  showing  the  two  Gaussian-like  components  and 
their  sum  from  the  measurements  of  Ref.  32.  The  maximum  value  is  set  to  1.0. 

Figure  A2. - fo(v),  O-atom  speed  distribution  distribution,  Eq  Al,  times  the  speed  squared, 

where  the  speed,  v,  for  fo(v)  refers  to  the  O-atom  speed  in  the  NO2  the  center-of-mass  frame; 

- fo-HCi(v),  relative  O-HCl  speed  distribution  from  Eq.  A2,  where  the  speed,  v,  for  fo-HCi(v) 

refers  to  the  O-atom  relative  velocity  in  the  O+HCl  center-of-mass  frame.  Note  that  each 
distribution  is  nonnalized  to  unity. 

Figure  A3.  O+HCl  relative  (normalized)  translational  energy  distribution  corresponding  to  the 
speed  distribution,  fo-HCi(v),  shown  in  Figure  A2. 
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FIGURE  9 
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FIGURE  10 
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FIGURE  A-l 
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FIGURE  A-3 
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